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ABSTRACT 

Aims. Rotation periods of cometary nuclei are scarce, though important when studying the nature and origin of these objects. Our aim 
is to derive a rotation period for the nucleus of comet C/2004 Q2 (Machholz). 

Methods. C/2004 Q2 (Machholz) was monitored using the Merope CCD camera on the Mercator telescope at La Palma, Spain, 
in January 2005, during its closest approach to Earth, implying a high spatial resolution (50 km per pixel). One hundred seventy 
images were recorded in three different photometric broadband filters, two blue ones (Geneva U and B) and one red (Cousins I). 
Magnitudes for the comet's optocentre were derived with very small apertures to isolate the contribution of the nucleus to the bright 
coma, including correction for the seeing. Our CCD photometry also permitted us to study the coma profile of the inner coma in the 
different bands. 

Results. A rotation period for the nucleus of P = 9. 1±0.2 h was derived. The period is on the short side compared to published periods 
of other comets, but still shorter periods are known. Nevertheless, comparing our results with images obtained in the narrowband CN 
filter, the possibility that our method sampled P/2 instead of P cannot be excluded. Coma profiles are also presented, and a terminal 
ejection velocity of the grains v gr = 1609 ± 48 m s is found from the continuum profile in the I band. 

Key words. Comets: individual: C/2004 Q2 (Machholz) 



1. Introduction 

Comets are at the same time both the largest and the smallest ce- 
lestial bodies in our solar system. While the nucleus of a comet 
is typically 1 to 15 km, the ion tai l can easily extend over more 
than 1 AU (e.g. iJones et al.ll2000h . Cometary nuclei are chem- 
ically and morphologically complex bodies, consisting of non- 
volatile ices and frozen grains. When these nuclei approach the 
sun, gases start to sublimate from them, forming the cloud of 
gas and dust known as the coma. For comets in the inner solar 
system, this coma is much brighter than the nucleus, making the 
nucleus hidden from our view. Yet, information on the dynamics 
of this nucleus is important to study not only the nature and ori- 
gin of these objects, but also more generally the origins of our 
solar system itself. In the first place, it is of a statistical impor- 
tance to determine cometary periods, since a representative sam- 
ple of periods could give important clues on possible statistical 
dependencies with other properties, like orbital period or mass, 
and with similar objects like asteroids. Moreover, information on 
the rotational state of a comet gives important clues on the recent 
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Table 1. Basic properties of comet C/2004 Q2 (Machholz). 



Comet C/2004 02 (Machholz) 



Discovery 


August 27, 2004 


Earth passage 


A = 0.35 AU on January 5, 2005 


Perihelion 


r h = 1.21 AU on January 25, 2005 


Orbital period 


-117000 year 


Inclination 


38.6° with ecliptic 


HtO prod, rate 


2.6 10 2S molec. s" 1 at A = 0.394 AU (t) 



Refs.: Horizons On-Line Ephemeris System of JPL 



(NASA) except f: Bonev et al. 2006 



history of the comet; fast rotation for example can point to strong 
recent activity in which the comet changed its own rotation state, 
or alternatively (but less likely) to a recent interaction with other 
solar system objects. Finally, rotational periods are also an im- 
portant jwarnejterjnjhe development of sublimation models for 
nuclei (ISchleicher et al.lll998l) . 

In this paper, we present a study of the Oort Cloud comet 
C/2004 Q2 (Machholz). Some basic properties of this comet are 
summarised in Table [ T] It was discovered o n August 27, 2004 
by Don E. Machholz dMachholz et al.ll2004l) . and reached a vi- 
sual brightness of ~3.5 mag. at its closest approach to Earth, 
making it one of the most spectacular astronomical events of 
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that year. The observations presented in this paper are taken 
around this close encounter, providing a high spatial resolution. 
In this paper, we concentrate both on the nucleus, and on the in- 
ner (10 5 km) coma. Our observational data consist of CCD im- 
ages in three photometric bands. Our first objective is to derive 
a rotation period of the nucleus from our intensive photometric 
monitoring, while our second goal is to study the profiles of the 
inner coma. These profiles are relatively simple to determine, 
and directly lead to some interesting properties, like the terminal 
ejection velocity of the dust grains. 

Published rotation pe riods of cometary nuclei are scarce (see 
ISamarasinha et al.1l2004 for a summary), since the nucleus of a 
comet situated in the inner solar system is always hidden in an 
opaque and overwhelming coma. Several techniques have been 
developed, each with their own drawbacks. One possible tech- 
nique is to search for structures in the inner coma that emerge on 
the nucleus, like fans or jets, and search for periodicities in their 
morphology. This technique was employ ed to search for the ro- 
tation state of C/2004 Q2 (Machholz) by iFarnham et al.1 d2007l) 
and will be discussed later in this paper (Sect. |5.3l l. The tech- 
nique used in this paper is aperture photometry using extremely 
small aperture radii. The idea is that in the vicinity of the op- 
tocentre of the coma, a small but observable contribution of the 
nucleus to the total b rightness is expected. In the formalism of 
iLicandro et al.1 J2000h . if the brightness measured in an aperture 
p around the optocentre is given by 

B(p) = B N + B c (p) 

where Bcip) is the brightness of the coma captured in a radius p 
and Bn the nucleus brightness. The amplitude of the light curve 
is then 

Am(p) = -2.5 log[min(fl w ) + B c (p)]/[max(B N ) + B c (p)]. (1) 

Note that in this equation, we implicitly assume that Bcip) is 
constant during the time span of the observations. In practice, 
changes in the activity can imply large coma variabilities, es- 
pecially near perihelion. Nevertheless, from Eq. dT|i, it is easily 
seen that when taking large apertures, the contribution of the 
coma Bcip) cancels out the amplitude of the light curve induced 
by the variability of the nucleus. Therefore, very small apertures 
around the comet's optocentre ar e usually employed to detect 
variability induced by the nucl eus dJewitt & Meechfl985t iJewittl 
ll990HMeech et al.ll993lfl997t) . However, a problem arises when 
taking apertures smaller than the seeing disc of the obse r vation . 
This problem was studied in detail by ILicandro et al.l (120001) . 
who dubbed this problem the seeing effect. They found that sev- 
eral published periods were in fact contaminated, or even at- 
tributable, to this seeing effect. In our study, we will take this 
effect into account. 

The reduction of the data, as well as its correction for the 
seeing effect, are an important part of our work, since it will 
ultimately determine the quality of our analysis. Therefore, the 
technical part of our analysis will be a considerable fraction of 
this paper: The observations are described in Sect. [2] while the 
reduction is discussed in Sect. |3j Sect. |4]is devoted to the extrac- 
tion of the cometary magnitudes, including the correction for the 
seeing effect; Next, the variability analysis and its results are pre- 
sented in Sect.[5j In Sect. [6] we construct the coma profiles in the 
three filters that we used for this programme. The general dis- 
cussion of the obtained results is found in Sect. |7J We end with 
the conclusions (Sect. [8}. 



Table 2. Observational log. For each night, the number of obser- 
vations is given for each filter U, B and I. Also the heliocentric 
(r;,) and geocentric (A) distances are given (in AU, at 0:00 UT), 
together with the correction factor/,, (see Eqn.|2]in Sect. |4.2| >. 

Observational log 





n(U) 


n(B) 


n(I) 


ri, 


A 


fa 


2005-Jan-08 


5 


5 


5 


1.2337 


0.3485 


1.00 


2005-Jan-09 


9 


9 


9 


1.2305 


0.3500 


1.00 


2005 -Jan- 10 


10 


9 


8 


1.2274 


0.3520 


1.01 


2005-Jan-ll 


8 


9 


8 


1.2245 


0.3546 


1.02 


2005-Jan-12 


8 


8 


8 


1.2219 


0.3577 


1.03 


2005-Jan-13 


8 


8 


8 


1.2194 


0.3613 


1.05 


2005-Jan-14 


6 


6 


6 


1.2171 


0.3654 


1.07 


2005-Jan-15 


5 


5 


5 


1.2150 


0.3699 


1.09 



2. Observations 

Comet C/2004 Q2 (Machholz) was intensively monitored with 
the 1.2 m Mercator telescope at Roque de los Muchachos at 
La Palma (Spain), between January 7 and January 15, 2005. 
The telescope is equip ped with a photometric camera Merope 
dDavignon et al.ll2004 . consisting of a 2kx2k EEV CCD with 
a filt er wheel consisting of the traditional 7-band Geneva sys- 
tem dGolavlll966l;lRufener & Nicoletl ll988l). added with a Gunn 
R dThuan & Gunnll 19761) and a Cousins I dCousinslll975l) filter. 
The field of view is a square 6/5, resulting in a 0."19 sampling. 

181 Images were taken, in three different bands: 61 observa- 
tions in both the Geneva U and B band, and 59 in the Cousins I 
band. A log of the observations is given in Table |2j while three 
typical observations are shown in Fig. [2] Integration times were 
fixed at 60 sec for B and I exposures, and 600 sec for U band ex- 
posures. Each set of observations was accompanied with a stan- 
dard star observation, immediately before, after or inbetween the 
comet observations. Standar d stars were tak en from the Geneva 
catalogue of standard stars dRufenerlll999l) . Special care was 
taken to ensure that the standard stars were not situated in the 
extended coma or tail of the comet. Standard stars are used to 
accurately determine the extinction conditions of the time of ob- 
servation. In addition, 40 sky flats in the three used filters were 
taken in this observing run, 16 in U, 11 in B and 13 in I. Few 
comet observations (~10) could not be used for further analy- 
sis, mostly due to a wrong tracking or pointing of the telescope. 
Weather conditions were exceptionally good with clear nights 
and a superb seeing throughout the whole observing run. 

Before we continue in discussing the actual reduction of the 
CCD frames, it is a very useful exercise to get an idea which 
cometary emission these three filters include. For this purpose, 
we downloaded low-resolution optical spectra of C/2004 Q2 
(Machholz) from the ESO archive (http://archive.eso.org). 
These archive data are taken on December 9, 2004, with the 
EFOSC2 instrument that is mounted on the ESO 3.6m tele- 
scope, under observation programme 274.C-5020 (RI. Weiler). 
Reduction was done in the available ESO-MIDAS routines by 
ourselves, and included the standard steps of a long-slit re- 
duction. Extraction was done over the whole spatial direc- 
tion. We succeeded to identify all major emission features 
in these spectra, taking laboratory wavelengths from different 
sources dKawakita & Watanabd 120021: lLara et all l200l lo^deil 
19711). High resolution atl ases of cometary emission lines (e.g. 
Cochran & Cochran] 120021) indicate that many more emission 
features are present, but these are not resolved at this relatively 
low resolution. From a simple inspection of this figure, it is clear 
that the Geneva U and B bands are sampling both continuum ra- 
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Fig. 1. The low-resolution spectrum of comet C/2004 Q2 (Machholz) taken with the EFOSC2 instrument taken with the ESO 3.6 m 
at La Silla on 9 december 2004. See the text for the line identification. The passbands of the three photometric filters of the Merope 
instrument used in this paper are shown in red. It is obvious that these broad filters cover a mix of continuum and emission features, 
with the I filter having the most emission-free passband. The two spectra are ESO archive data, and are taken under programme 
274.C-5020 (P.I. Weiler). 




Fig. 2. Three observations, one in each band, taken on January 
11, 2005. The frames are shown here to give an idea of the type 
of data. They are taken consecutively, and hence under very sim- 
ilar ambient conditions. In the U filter, the very clear signature 
of the ion tail is visible, in B this tail is only marginally visi- 
ble, and in I it is absent. The field of view of these images is 
3/6 x 5/0 of the total square 6/5 of the CCD; the scale is shown 
in the lower left corner. Note that on these images, the core of 
the coma seems to be saturated, which is obviously not the case, 
but solely a consequence of the chosen colour scale. 



diation and emission features: NH (and most probably also the 
very strong OH emission at 3085 A) in the U band; in the B band 



the strong CN (Av = 0) feature, but also C3 and CH features. The 
I filter, however, does not contain such strong emission features, 
and mainly samples the dust continuum. 

3. Reduction 

The reduction of CCD frames involves the bias correction, the 
construction of a master flat field, the correction of the raw sci- 
ence frames by this master flat field, and the correction for the 
background. 

Bias correction was done by using the pre- and overscan re- 
gions of the image. It was not possible to construct one master 
flat per filter for the entire observing run, since a few dust par- 
ticles situated in front of the CCD were clearly visible on the 
flats. Therefore, several master flats were created per filter, each 
of them constructed with raw flats that show the same pattern 
of dust particles. Flats were combined by a median and sigma- 
clipping method, while special attention was drawn to an ade- 
quate removal of the few weak stars present in some raw flats. 
Note that none of the pixels that were affected by the removal 
of those weak stars, coincide with pixels that were used in the 
extraction of the comet's optocentre (Sect. |4~4t . 

The background correction was also not straightforward. In 
the case of the standard stars, a clear gradient is present in the 
background, probably caused by the electronics involved in the 
windowed (300x300 pixels) read-out of these frames. The gra- 
dient was, however, efficiently removed by the construction and 
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Fig. 3. Correction for the seeing effect for a frame taken in opti- 
mal conditions (left column of figures), and for a frame taken 
in less favourable conditions (right column). The upper row 
shows the optocentres of the reduced frames (but not yet seeing- 
corrected). To degrade these frames to the same seeing, the left 
one has to be convolved with a broad Gaussian kernel (middle 
left panel), while the right one with a very narrow Gaussian ker- 
nel (middle right panel). The resulting optocentres are show in 
the bottom panels. 



subtraction of a second degree polynomial. In the case of the 
comet frames, the background could not be determined on the 
frame, since the coma of the comet is much larger than the field 
of view of our detector. No dedicated sky frames were taken dur- 
ing the observing run. Consequently, no background correction 
was done on the comet frames. This choice can be motivated by 
two arguments. First, a rough estimate of the sky contribution 
could be made on science frames of other science programmes 
that night, and we concluded that on all frames the sky contribu- 
tion is several orders of magnitude weaker than the coma, even 
on the borders of the CCD. Second, if we fit the decline of the 
coma with a polynomial f(p) of nth degree, starting from the 
centre of the coma, then the limit of f(p) for p— >oo is a useable 
estim ate of the background contribution (e.g. lSchulz & A'Hearnl 
fl995l) . In all our fits, this limit do is significantly small, al- 
though there is some difference in «o for different polynomial 
degrees, and for different fitting directions away from the centre. 
Moreover, the observations were carried out around new moon 
(January 10), assuring a minimal sky contribution. In any case, 
the background does certainly not exceed a few thousandths of 
the pixel values around the centre of the coma. Therefore, back- 
ground subtraction can be safely neglected when concentrating 
on the central luminosity, as we do in this paper. 



4. Extraction 

4. 1 . Extraction of standard stars 

Due to their brightness, standard stars are observed out 
of focus. Standard stars are extracted using SourceExtractor 
dBertin & Arnoutsl I1996I) . The method of extraction is a sim- 
ple addition of the flux above a certain threshold (1.4cr), with 
a Gaussian correction for the wings of the star. From these stan- 
dard star magnitudes, extinction coefficients kx are derived in 
each filter. These extinction values were then interpolated to the 
observation times of the comet. 

The extinction values in the U and B filters correlate very 
well, which is an independent check for their accuracy. The 
mean kx, the range (1 cr) and the peak-to-peak of the extinction 
values in the U filter for the whole observing period are as fol- 
lows: kx = 0.55 ± 0.03 with peak-to-peak 0.12. For the B filter, 
these values are kx =0.30 + 0.02 and peak-to-peak 0.07. These 
basic statistics show that the extinction values are very stable 
throughout the observing run. We have also performed a fre- 
quency analysis of the extinction values for these two filters, but 
no significant frequencies were found in the range [0, 10] cd . 

The extinction values in the I filter are much more variable 
than the values for U and B. Moreover, they do not correlate with 
the two other filters. Since we could not find an explanation for 
this result, we chose to derive the extinction in the I band from 
a linear combination of the U and B band extinctions. Luckily, 
the extinction correction is not so critical for the I band, since it 
is generally very small ($0.1 mag). 

4.2. Distance correction of the comet frames 

Before the actual extraction of the coma centre, one has to cor- 
rect for the motion of the comet both relative to the observer and 
the sun. During our observing run, both effects have an opposite 
effect, as the comet approaches the sun while it recedes from the 
earth. To correct for the luminosity change caused by the motion 
of the comet, we have multiplied the science frames by 



with A the distance comet-earth, and r;, the distance comet- 
sun, whereas the subscript denotes these distances on the 
first night of our observations. The values of these quan- 
tities were retrieved from the Horizons On-Line Ephemeris 
System from the Solar System Dynamics Group of JPL (NASA) 
(http ://ssd . jpl . nasa . gov/horizons . cgi] > and are given in 
Table [2] The factor f a reaches a value of 1 .09 in the last night, 
resulting in a magnitude decrease of ~0.04. 

4.3. Correction for the seeing effect 

In this particular case, especially the extraction of the science 
frames (see Sect. 14.41 ) is a delicate step, since we will deal 
with extremely small aperture sizes, much smaller than the see- 
ing disc. Since seeing conditions are expected to be variable 
throughout the night, extraction with a aperture size smaller than 
a typical seeing will introd uce artificial va r iabilit y. In a paper de- 
voted to this seeing effecf . ILTcandro etaD (120 00*) claim that sev- 
eral published rotational periods of comets are doubtful, and are 
probably affected by this seeing effect. It is therefore extremely 
important to account for the variation in the seeing. 

Detailed and accurate seeing measurements are available on 
the Roque de los Muchachos site, provided by the dedicated 
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RoboDIMM telescope dAugusteiinll200U lO'MahonvlfeOCBl) of 
the ING. RoboDIMM forms four images of the same star, mea- 
suring image motion in two orthogonal directions from each of 
the two pairs of images, from which it derives four simultaneous 
and independent estimates of the seeing. The values are auto- 
matically corrected for the zenith distance and a wavelength of 
550nm. The time resolution of these measurements is 2 minutes, 
and therefore very well suited for our project. We will use these 
data for an adequate correction of the seeing effect. The idea is to 
degrade the seeing of each science frame to the worst seeing that 
was recorded during the comet observations. While this method 
seems to be quite rough, its application can be motivated by two 
facts. First, thanks to the good seeing conditions, the "worst see- 
ing" was not that bad at all. In Fig |4] the RoboDIMM seeing 
values interpolated to the comet observation times are shown. It 
is clear that the least favourable seeing does not exceed a FWHM 
of ~2."2. Second, the reverse process, in which frames taken in 
bad seeing conditions are deconvolved to an average or a very 
good seeing, is a much less robust process. The comets frames 
were degraded to the worst seeing in that filter, by convolution 
of a twodimensional Gaussian kernel. As a consequence, images 
taken in good seeing conditions are convolved with a broad ker- 
nel, while images taken in worse conditions, only need a small 
kernel to be degraded with. On Fig. [3] two examples are show 
of such a convolution. The left column shows the technique for 
a frame taken in good conditions, while the right column shows 
the same set of frames for a frame taken in bad conditions. From 
top to bottom this figure shows: the frames as they come out 
from the reduction (upper row), the Gaussian kernel to convolve 
with (middle row), and the resulting, convolved frames (lower 
row). The upper left frame shows the optocentre of a comet 
frame taken at good conditions. Consequently it needs quite a 
broad kernel (middle left frame) to degrade it to desired seeing. 
The reverse applies to the right set of frames, with the raw frame 
taken in less good conditions. 

Note that an objection could be made against the use of 
seeing data originating from another source than our own tele- 
scope. Indeed, the seeing is computed for a region of the sky 
that is not necessary the same as the region containing the comet. 
Nevertheless, we chose this option for several reasons, (i) A large 
part of the observing run was photometric, and hence only small 
spatial variations of the seeing value are expected; (ii) Our own 
observational data do not permit to deduce a seeing value: the 
standard stars are taken out of focus (to prevent saturation), the 
science frames are tracked on the comet, and no other science 
frames are available inbetween the comet observations. 



4.4. Extraction of the comet 

Now that all comet frames are degraded to the same seeing, we 
can safely extract the optocentre of the comet with small aper- 
tures well below the seeing disc. This technique poses, on its 
own turn, two specific problems to deal with. The first one is 
the determination of the exact location of the optocentre, while 
the second one involves the discretisation of the pixels when us- 
ing very small apertures. The extraction is then done using three 
different extraction radii (1,3 and 6 pixels). 



4.4.1 . Determination of the centre 

Several solutions can be put forward to determine the optocentre 
of the coma. The most simple one is obviously taking the po- 
sition of the brightest pixel in the coma. A second one is to fit 



10 12 
Julian Date (- 2453372) 



Fig. 4. The evolution of the seeing during ou r observational run 
as recorded by the RoboDIMM telescope (lAugusteijnl l200lb 
lO'MahomteCKBl) . interpolated to the times of our comet obser- 



vations. 



a twodimensional Gaussian in the proximity of the optocentre. 
A third possibility is to take the barycentre of a limited region 
around the optocentre. We chose this last method, since some 
tests showed that the first method was a too rough determina- 
tion, while the Gaussian method did not give good results either. 
The Gaussian method obviously assumes a Gaussian profile for 
the optocentre. While this is indeed the case for the immediate 
surroundings of the optocentre due to the effect of the seeing, it 
is not the case for a larger area since the coma shows a profile 
different from Gaussian (see Sect. [6j. 

4.4.2. Discretisation of pixels 

In the case of extraction with small apertures, flux has to be 
added from the sub-pixel level. Common extraction programs, 
such as SExtractor, divide each pixel in a certain number of vir- 
tual sub-pixels (5x5 in the case of SExtractor), and add the flux 
in those sub-pixels lying in the predefined aperture radius. Here 
we chose to work strictly analytic. For each pixel lying at the 
border of the aperture radius, we calculated the exact fraction of 
that pixel in the radius. 



5. Variability analysis 

The final magnitudes for the optocentre of the comet are 
graphically presented in Fig. [5] A detailed frequency analy- 
sis was performed on this data set, using tw o different algo- 
rithms: the Jurkevi ch-Stellingwerf pdm met hod ( Jurkevichl 1971 



Stell ingwerfll978h . and the Lomb-Scargle (lLomblll976HScarglel 
1983) periodo gram. The latter one is essentially a Fourier tech- 
nique, hence it is more sensitive for harmonic periodicities. The 
former technique searches for the frequency for which the peri- 
odogram, devided into a certain number of bins, shows the small- 
est spread per bin. 

As a first step, a frequency interval has to be defined. A 
conservative choice would be the interval between zero and 
the Nyquist frequency. In our case, when taking the inverse of 
the average time gap between neighboring points while ignor- 
ing large gaps, the value of the Nyquist frequency is 44cd~' 
(P = 0.55h). Such a high frequency will, however, lead to a 
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Fig. 5. Overview of the magnitudes for the extracted optocentre 
of the comet in the three different bands. Shown here are the 
results for an aperture radius of 1 and 6 pixels (0." 19 and I." 14). 



break-up of the comet. A rough estimate of the break-up fre- 
quency fb can easily be calculated by stating that at the equa- 
torial velocity vj the centrifugal acceleration equals the grav- 
itational acceleration, and that there is no internal strength to 
the body. If one assumes a typical density of lgcirT 3 , then 
the break-up frequency will be fb = 7 .3 cd~ l (P/, = 3.3h). For 
smaller densities, the break-up frequency wi ll obviously be 
small er; the same is true for aspherical nuclei dJewitt & Meechl 
1 19881) . Here, we adopt a safe upper limit for our frequency inter- 
val of / = 20 c d ~' (Pb = 1.2 h). The cometary rotation f requen- 
cies that are found up to now dSamarasinha et al]|2004l) all fall 
safely in this range. 

5.1. Non-periodic variability 

As already seen from Fig. [5] there is a remarkable increase in 
brightness in the nights 12-13 and 13-14 January of the observ- 
ing run. The brightness increase in those two nights leads to a 
frequency of / = 0.20cd -1 (P = 5d) in the Lomb-Scargle peri- 
odograms (Fig. [6), in each filter, and for each of the three radii (1, 
3 and 6 pixels). This frequency is also present in the Jurkevich- 
Stellingwerf results, but somewhat less clear. 

The increase in brightness is also seen in the "raw" magni- 
tudes, i.e. the magnitudes that are not seeing-corrected. Hence, 
the brightness increase is certainly not an artifact of our convolu- 
tion technique explained in Sect. 14.21 However, since the convo- 
lution technique is certainly not perfect, it might still be possible 



that the brightness increase is related to the seeing values that are 
somewhat higher than average for these two nights. We studied 
this possibilty by making a variability analysis of the seeing val- 
ues, and concluded that none of the variability that was found in 
the seeing, was still present in the seeing-corrected magnitudes, 
while the raw magnitudes showed the same trends as the seeing 
values, as expected. Hence we are confident that the seeing effect 
is adequately accounted for. 

Since the total time span of our observational run is only 
8 nights, it is not possible to judge whether this frequency is 
real (implying that the brightness increase is a periodic phe- 
nomenon), or it is just a unique event. In any case, since 
the typical rotational frequencies of cometary nuclei are much 
higher, it is unlikely that this frequency is linked to the rota- 
tion of the nucleus. The brightness increase is more likely a 
non-periodic event, like an outburst. Possibly it is connected to 
the s udden change in the solar wind velocity around that pe- 
riod dDegroote et alj|2008l) . It is, however, unlikely that the solar 
wind has a direct impact on the comet's nucleus, since it is be- 
lieved that s olar wi nd cannot enter the so-called contact surface 
(<10 3 km,||3|2004|) of the comet. 



5.2. Rotation frequency 

In order to search for an acceptable rotation frequency, the 
brightness increase should be removed. This can be done in 
several ways: e.g. by means of removing ("prewhitening") the 
data with the 5 days periodic component found earlier. Another 
possibility is to subtract the mean value of every night. Since 
these means are often ill-defined because of the scarcity of the 
data, we tried to model the data with a combination of unknown 
nightly means and a periodic component. A third possibility is 
to perform a simultaneous fit of two frequencies together, one 
long-term component and one short-term. We will briefly dis- 
cuss these three methods and their results. 



5.2.1. Prewithening 

When the frequency / = 0.20 c d -1 (P = 5 d) was removed from 
the data, we did not find any significant frequencies anymore 
in the B filter magnitudes. Two different frequencies pop up in 
the U and I filter periodograms. In the U filter, the most sig- 
nificant frequency is present near / = 6.03 c d -1 (P =4.0h), but 
after inspection of the periodogram, this frequency and some 
other frequency close to integer values, seem to be aliases of the 
frequency led -1 . There is little doubt that this frequency is an 
artifact, probably introduced by the imperfect extinction correc- 
tions. Moreover, the probable contribution of OH emission in the 
U filter (see Sect. |2]i could explain why this one-day alias is so 
prominent for this filter. Indeed, the extinction at OH is very high 
(~2mag/airmass), so as the comet rises or sets, the signal from 
OH can change dramatically. In the I filter, where overall extinc- 
tion is much less, a significant frequency is seen at / = 2.64 c d -1 
(P = 9. 1 h) in both the 1 and 3 pixel size apertures; for the 6 pixel 
magnitudes, a frequency of / = 3.64cd~' (P = 6.6h), which is 
a one-day alias of the former one. The frequency of 2.64 cd l 
(P = 9.1h) is also found in this filter and apertures when using 
the Jurkevich-Stellingwerf method. The variance reduction by 
this frequency is another 14% (from 50% to 64%). 
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5.2.2. Free nightly means 

After subtracting the mean values of every night from the 
data, the results of the frequency analysis around 2.5 cd~' (P 
around 9.5 h) were not different from the results obtained after 
prewhitening. However, some spurious, mostly higher, frequen- 
cies are present as well. They appear because the nightly mean 
values are not necessarily the true mean values of the light curve, 
since the rotation period could be longer than the time span of the 
comet observations in one night and since only few observations 
per night were available. A solution to this is to perform a har- 
monic fit through the data, while taking the mean value for each 
night as a free parameter. This technique gave no results for the 
U filter data, but now a frequency near / = 2.50 c d _I (P = 9.6 h) 
appears in the B filter data and near / = 2.59 c d -1 (P = 9.3 h) in 
the I filter. These results should be looked at with some caution, 
since 1 1 unknowns (three for the harmonic fit, and one for ev- 
ery night) were to be determined in data sets of only about 50 
observations. 

5.2.3. Two frequencies simultaneously 

Finally we looked also for the best two-frequency fit in the data. 
We searched for the best description of the data in the least 
squares sense with one long-period component (/ between 0. 1 
and 0.3 cd~', corresponding to P between 3.3d and 10 d) and 
a shorter-period component (/ between 0.3 and 3.0cd _1 , cor- 
responding to P between 8h and 3.3 d). The f — 1.00c d -1 fre- 
quency and its aliases are still too dominant in the U filter data to 
give reliable results. In the I filter data, a long term component of 
0. 18 c d -1 (P = 5.6 d) combined with the frequency / = 2.67 c d _I 
(P = 9.0h) reduces the variance of the data with 65%. In the 
B data the combination of 0.22 cd" 1 (P = 4.5d) and 2.56 cd~' 
(P = 9.4h) describes the observations best and leads to a vari- 
ance reduction of 83%. 

5.2.4. Final rotation frequency 

A summary of the frequency analysis is given in Table [3] In this 
table, the results are given in cycles per day (cd _1 ), since this 
notation allows us to immediately detect one-day aliases, and to 
compare directly with the periodograms of Fig. [6] The results 
of the different methods indicate that a short-term periodicity is 
present in the data, with a period between 9.0 and 9.4 h. Since not 
all methods gave consistent results for the B filter data, probably 
also due to imperfect extinction correction, we take the value 
found in the I filter data (9.1 h on average) as the best estimate 
for the rotation period from these data. The spread of the results 
given by different techniques, indicates that the error (1 cr) on 
this value can be as high as 0.25 h. 

5.3. Rotation of Machholz in literature 

ISastri et al.l (120051) analysed R band images of C/2004 
Q2 (Machholz) taken in the same period as our observa- 
tions (January 2005), and report the same period as ours: 
P = 9.12±1.9h (/ = 2.63cd _1 ). The period was derived by 
studying dust fans visib le in these sm all band images, and adop- 
tion of Fulle's formula dFulldl 19871) between grain velocity and 
the forces on each grain. It is a remarkable result that two such 
different methods lead to the same value for the rotation period, 
str engthening of course th e confidence in our results. 

IFarnham et all d2007l) observed C/2004 Q2 (Machholz) in- 
tensively during several observing campaigns from January to 



Table 3. Summary of the frequency analysis results. Three dif- 
ferent methods were used to get rid of the brightness increase 
during two nights of our observational run (12-13 and 13-14 
January); see Sect. l5.2l for more details. 



Frequency analysis results 


#1: 


after prewithening with f = 0.20 c cT l 


U 


strong one day aliases 


B 


no significant / 


I 


2.64c d" 1 (1, 3pix) 




3.64 cd~' (6pix) 


#2: 


after subtracting (free) nightly means 


U 


no results 


B 


2.50 cd" 1 (1, 3pix) 




0.50 cd" 1 (6pix) 


I 


2.59 cd _1 (allpix) 


#3: two frequency fit (simultaneous ) 


U 


strong one day aliases 


B 


0.22 c d" 1 and 2.56 c d" 1 (all pix) 


I 


0.18 cd" 1 and 2.67 cd" 1 (allpix) 




frequency (c/d) frequency (c/d) 



Fig. 6. Jurkevich-Stellingwerf pdm (left panels) ^-statistic and 
Lomb-Scargle (right panels) periodogram for the I magnitudes 
of 1 pixel aperture. In the two upper panels, the diagrams for the 
original data are shown, while in the lower panels, the diagrams 
for the magnitudes prewhitened with f\ = 0.20 c d -1 (P = 5 d) are 
shown (method #1 in Table [3). On each figure, the frequency 
/ = 2.64 c d" 1 (P = 9. 1 h) is shown with a red line. 

June 2005, in different filters of the HB narrowband comet fil- 
ter system (IFarnham et al.ll2000h . The paper concentrates on the 
CN band images, which reveal two clear jets in roughly opposite 
directions relative to the nucleus. The presence of the jets were 
confirmed by images taken with the 1-m telescope at the Lulin 
Observatory, Taiwan, and th e 3.6m at ESO, La Silla dLin et al.l 
120071) . IFarnham et alJd2007l) found that the morphology of these 
jets repeated itself on the timescale of hours, and through an in- 
tensive monitoring in April 2005, the authors could phase the jet 
morphology on a period of 17.6±0.05 h, which is then also iden- 
tified as the rotation period of the nucleus. fFarnham et al.l d2007l) 
explicitly excludes the possibility of a period shorter than 17.6 h, 
since they did not find similar morphologies of the CN jets with 
a time span shorter than this period. Also their March obser- 
vations are consistent with this 17.6 h period. They argue that 
it is unlikely that the rotation period of the nucleus should have 
evolved from 9. 1 h to 17.6 h between January and March. Hence, 
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Fig. 7. Phase diagrams of the comet's optocentre magnitudes in 
the I filter with an aperture radius of 1 pixel, upper panel: phase 
plot with f\ =0.20cd~' (P = 5d); lower panel: phase plot with 
/2=2.64c d -1 (P = 9.1 h) on the magnitudes prewhitened with f\ 
(method #7 in Table 0. 



the period of ISastri et all as well as the peri od found by us , is 
seemingly in contradiction with the one from lFarnham et aU 

There might be, however, a possible explanation th at could 
reconcile the two results. The method of iFarnham et al.l is based 
on the visual inspection of the morphology of the CN jets that 
emerge from the nucleus. Since it concerns two opposite jets, 
there is always one jet dominant in the direction of the ob- 
server. If we assume the jets, or more precise, the regions on 
the comet's nucleus where the jets emerge, to be the main con- 
tributors in the comet's variability, then it is obvious that our 
method c atches P/2 and n ot P. If this interpretation is correct, the 
re sults of Farnham e t al.l are in agreement with ours and the one 
of ISastri et al.l Note that the latter period was derived by mod- 
elling the structure and position of the observed dust fans, and 
not by a period icity search in a series of images, as was done by 
IFarnham et al.L 

In the assumption that our method leads to an estimate of 
P/2 and not P, even then a sma ll difference persists between our 
P/2 = 9.1 h and IFarnham et all 's P/2 = 8.8 h. This difference can 
have different explanations. Surely, the measurement errors will 
account for an important part in this difference, but also another 
effect can be of importance. Table [2] shows a 7% change in A 
from the first to the last night's data. This produces a commen- 
surate change in the aperture sizes when measured in km at the 
comet. If the variability is caused by changes in the comet's ac- 



tivity, then the period determination could be affected by this 
difference. For example, a burst of activity is seen as a bright- 
ening of the coma. When the leading edge of that burst leaves 
the aperture, the light curve has reached its peak and starts fad- 
ing again. Because a given aperture is smaller (in km) on earlier 
nights, the observed peak will occur earlier than it should rela- 
tive to the last night. This produces a period measurement that 
is longer than the actual periodicity (defined by the burst of ac- 
tivity reaching the same distance from the nucleus). Given the 
small pixel dimensions at the comet and the high emission ve- 
locities, this could have a significant effect, an d thus may hel p 
to explain the difference between our 9.1 h and IFarnham et all s 
8.8 h half -periods. 

6. Coma profiles 

The first objective of our observational programme was to study 
the rotation of the nucleus through extremely small aperture 
photometry. The use of the CCD camera implies, however, that 
we sample a large part of the inner coma at each observation. 
Therefore, we decided to initiate a study of the coma profile in 
the three spectral bands at our disposal. In Sect. |6.2| we will de- 
duce some physical properties of the coma from these profiles, 
but first a small introduction is given of what can be expected 
theoretically. 



6. 1 . Theoretical profiles 

Historically, the fi rst quantita t ive m odel of a come tary coma 
was elaborated by lEddingto"nl dl910h . In lEddingtonl 's fountain 
model the comet is assumed to be a uniform and isotropic source 
of emitting particles such that their density (e.g., gas or dust) 
would fall as the inverse square of the distance from the source, 
except the emitters are also subject to a uniform acceleration 
that pushes on the m from a given dir ection, pres umably from 
the S un's direction (Combi et al.ll2004l) . Continuing lEddingtorif s 
work. lWallace et aL ( 1958 ) showed that the column density N in 
the line of sight, and hence the brightness B in the case of an 
optically thin coma, is given by 



N = 6/(4vp) 

where Q is the global particle production rate, v is the initial 
uniform outflow speed and p is the projected distance on the 
sky plane of the nucleus. Generally the logarithmic derivative 
m — af logB(p)/iil ogp is used to describe the coma profile, and 
hence m--l in I Wallace et all s equation. Deviations from the 
m = - 1 profile are expected at distances where t he solar radiation 
pressu re becomes important. More specifically. iJewitt & Meechl 
(119871) showed that this distance is of the order of 



X R * v^/PAWl)] 



(3) 



where v gr is the terminal ejection velocity of a grain, g sun is 
the solar gravity at 1 AU, and B is the ratio of the acceleration 
of a grain due to radiation pressure to the local solar gravity. 
Beyond this poi n t Xr, one enters the outer coma where |m| > 1 . 
IJewitt & Meechl (1 19871) derive from both Monte Carlo simula- 
tions, as well as from an analytical derivation, that m^-l .5 in 
this region, slightly varying depending on the phase angle and 
the radiation pressure. 

The theoretical profiles that are discussed here assume that 
the coma has constant (though not necessarily isotropic) activity. 
In realit y, C/2004 Q2 (Machho lz) is known to have significant jet 
activity (IFarnham et al.l 120071) . and jets produce structures that 
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can affect the radial profiles in the coma. Jets on rotating nuclei 
frequently produce curved or spiral shaped jets. Radial profiles 
may cross these curved features, causing the inner or outer re- 
gion to be brighter than otherwise might be expected. This is es- 
pecially true in the inner coma we are measuring. Since activity 
changes are invoked as the cause of the variability that defines 
the cometary period (Sect. 15. 3K the profiles are expected to be 
influenced by this effect. 

6.2. Coma profiles of C/2004 Q2 (Machholz) 

To construct representative coma profiles of C/2004 Q2 
(Machholz), we took for each filter one observation taken in the 
night 08-09 January 2005, which was a night with optimal con- 
ditions and seeing (see Fig.|4|i. For each of these three observa- 
tions, we constructed two coma profiles, one in the direction of 
the sun, and one in the direction of the ion tail. The profiles were 
constructed by dividing the coma into concentric annuli, with 
radii that are (logarithmically) equidistant. For each of these an- 
nuli, two medians were taken: one of the sunward half of the an- 
nulus (sunward profiles in Fig. [8]) and one of the tailward half of 
the annulus (tailward profiles in Fig. [8]). The profiles were nor- 
malised so that the brightness of the optocentre reaches unity, 
and plotted on a logarithmic scale. The resulting profiles are 
shown in Fig. [8j and linear fits were made in the two regimes 
2.6 < log p < 4 and 4 < log p (p in km). The values for m are also 
given in this figure. 

While the profiles in the I are fully compatible with the ba- 
sic models presented in the previous section, the profiles in U 
and B are clearly deviating from these predictions. From Fig.Q] 
it is clear that in both the U and B bands strong emissions are 
present, while the I filter is relatively free from strong emission 
lines. The presence of these emissions breaks the assumption of 
conservation of mass of a certain species, since the gas can dis- 
sociate and/or ionise to other species. The development of theo- 
retical profile s taking these effects into account were initiated by 
iHaserl (119571) . and are specifically valid when using small band 
filters concentrated on those emission described in the profile. 
Since we are using broadband filters in this study, sampling both 
continuum and emission features, it is very difficult to construct 
a valid profile for the U and B filters. The I filter, however, is 
dominated by the continuum, leading to a clean m = - 1 profile 
in the inner coma. 

From the change in m from -1 to -1.5 at a distance of 
Xr =21.6+ 1.3 x 10 4 km, we can estimate the terminal ejec- 
tion velocity of the grains using formula ©. Typically, 
one assumes y8-l, which corresponds at 1 AU to spheri- 
cal grains with a diameter of about 1 pm. This way, we 
can estimate v gr = 1609 ± 48 m s -1 , which is much more than 
we expect from t he empirical law of Bobrovniko ff-Delsemme 
dPelsemmdl 19821) : v gl - = 580/ VPfe = 524 ms 1 (e.g. ISchulz etail 
Il993l) . Deviations from this law of this magnitude may be at- 
tributed to stages of high activity. For distances smaller than 
logp $ 2.6, there is also a clear deviation from the m = -1 trend. 
At this scale (p ^400 km), other effects come into play. The ac- 
celeration zo ne close to the nucleus h as certainly some influence 
at this scale dJewitt & Meechl[l987l) . but most importantly it is 
the seeing that wipes out the m = - 1 profile close to the optocen- 
tre. 

By calculating these ;7i-values for all observations, we can 
search for temporal variations or global trends in the profiles. 
We concentrate on the profiles in the I band because of their 
agreement with the theoretical profiles. We did not use observa- 
tions containing bright background stars for this purpose, since 



the profiles of such observations are heavily disturbed by the 
star tracks. For the determination of the inner and outer coma 
profiles, respectively 56 and 5 1 observations were used, and the 
fitted m-values are shown in Fig. [9] From these m-values, it is 
clear that the profiles show some variation indeed, but overall 
we can safely state that the profiles are quite stable throughout 
our observational run. The spread on the m-value for the inner 
coma for example, is m — —1.01 ± 0.03; for the outer coma this is 
m = -1.49 ± 0.08. The seemingly periodic behaviour of the inner 
coma profile (Fig. [9] upper panel) is very peculiar, as we cannot 
link it to another phenomenon: not with the global brightness 
trend we encountered previously, nor with the deduced period 
(P or P/2), nor with any solar wind events. The observed trend is 
also noticeable for the inner coma in the U and B band, but not 
in any of the m-values of the outer coma. The fact that we cannot 
match the variations in the slopes to any known periodicity might 
be explained if the slopes are being affected by coma structures, 
as discussed earlier. If the slopes are affected by temporary struc- 
tures, then sparsely sampling the slope at various times may not 
reveal variations that change with the rotation period. 

On Fig. [9] the evolution of the terminal ejection velocity of 
the grains v gr is also shown. This quantity is related to the helio- 
centric distance r/, of the comet, but the difference for r/, between 
the beginning and ending of our observational run is too small to 
see any noticeable effect. 

7. Discussion 

Rotation periods or periodicities are currently known for about 
25 comets. A comprehensive inventory of t he currently known 
or sus pected rotation periods is given in ISamarasinha et all 
d2004l). while a more c ritical list of 15 comets is given in 
iJorda & Gutierrez] (120001) . The mean of the spin periods given in 
IJorda & Gutierrezr s list is 20 h, but decreases to 12 h when con- 
strained to the short-period comets. Moreover, there seems to be 
a hint for a relation between rotation period and the size of the 
nucleus, with larger nuclei having longer periods. However, the 
small number of published periods does certainly not allow us to 
draw any final conclusions concerning statistical de pendencies. 

Th e difference between the period found by ISastri et all 
d2005l) that we conf irmed here, and the period fo u nd by 
IFarnham et all d2007l) is puzzling. ICrifo & Rodionovl d 19991) 
claim that the present circumnuclear models have too uncer- 
tain or even contradictory assumptions, implying that phys- 
ical properties for the nucleus derived from observations of 
jets or fans, can deviate considerably from reality. However, 
in the ca se of C/2004 Q2 (M achholz), the CN jets ob- 
served by IFarnham etalJ ( 120071) are bright, and the rotation 
is obvious from a simple inspection of the phased observa- 
tions (http://www.astro.umd.edu/~farnham/Machholz/). The 
authors reject a period shorther than 17.6 h, but, as explained in 
Sect. 15.31 there is a possibility that our method picked up P/2 in- 
stead of P. However, a somewhat uncomfortable inconsistency 
persists: if the periodicity of 9.1h that we found in the light 
curve, is induced by the periodicity in the morphology of the CN 
jets, then it is not clear why the period is best seen in the I band, 
and less clear in the B band (see Tabled, which contains the CN 
emission feature (Fig.[T). A possible explanation for why the pe- 
riodic variations are visible in the I filter, but not in the B filter, 
which captures the CN emission, is through the dust production. 
IFarnham et all d2007l) observed full CN spirals in C/2004 Q2 
(Machholz), indicating that they remain active throughout the 
rotation. Dust production, on the other hand, tends to shut down 
when a source passes into darkness. Thus, the brightness in the B 
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Fig. 8. Coma profiles of three representative observations for the 
three filters used in our programme. Crosses represent the profile 
in the sunward direction, while the xes show the profile in the 
direction of the ion tail. The coma in the I filter shows a clean 
m = - 1 profile for p < 10 4 km and m = -1 .5 for p > 10 4 km, the 
profiles in the U and B filter have smaller values for |m|, indicat- 
ing the influence of gas emissions in these spectral bands. 



Fig. 9. Temporal evolution (in function of the heliocentric dis- 
tance r/,) of the inner and outer coma profiles represented by 
their m-values, and of the terminal ejection velocity determined 
with formula (0). The terminal ejection velocity is related to 
the heliocentric distance r/, of the comet, but the difference of 
Ar/, ^0.02AU between the beginning and ending of our obser- 
vational run is too small to establish global trends. 



filter, which has a significant signal from CN, may remain fairly 
constant, while the I filter, which is d ominated by continuu m, re- 
flects the dust turning on and off. In lFarnham et alj d2007l) a fu- 
ture paper is announced in which also the observations of C/2004 
Q2 (Machholz) in the other HB narrowband comet filters are to 
be discussed. This could hopefully clarify the different results. 

It is clear that a light curve analysis of a comet's optocen- 
tre is hampered by several uncertainties. Particularly, the correc- 
tion for the seeing effect is a delicate step, and the introduced 
uncertainties are difficult to quantify. Our observational setup 
is, however, unique in the sense that our spatial resolution of 
~48kmpixek' is much higher than papers using the same tech- 
nique of aperture photometry. Therefore, it is not surpirising that 
the period is also recovered even when using a relatively large 
aperture of 6 pixels. 



8. Conclusion 

Comet C/2004 Q2 (Machholz) was monitored using the Merope 
CCD camera on the Mercator telescope at La Palma, Spain, in 
January 2005, during its closest approach to Earth. 170 Images 
were recorded in three different bands. Our main goal was to find 
a rotation period of the nucleus through aperture photometry us- 



ing very small apertures. This technique was already succesfully 
applied to other brigh t comets, but it was recently criticised by 
iLicandro et al.l d2000h . These authors studied the effect of aper- 
ture radii smaller than the seeing disc, and concluded that many 
periods published in literature that make use of this technique 
are contaminated by this seeing effect. In our analysis of comet 
C/2004 Q2 (Machholz), we took this effect into account by de- 
grading our images to the least favourable seeing. Luckily, since 
the ambient conditions during our observational run were very 
good, only the images with a superb seeing had to be degraded 
slightly. 

We searched in each filter and three aperture radii (1,3 and 6 
pixels) for frequencies between and 20 c d _1 (i.e. periods down 
to 1.2 h), and this with different algorithms. After taking into 
account that a long term trend i s in our data (which i s possibly 
linked to a solar wind event, see lDegroote et alj|2008[) . a period 
near 9. 1 ±0.2 h was found. While this rotation period is in perfect 
Sastri et al. ( 2005[). it clearly deviates 
Far nham et al. 



agreement with the one of 
from the 17.6 h found by 
methodology of both methods could explain the possibility that 
our method in fact finds P/2 and not P. 

Also coma profiles for the inner coma (log p < 4.6, p in km) 
were constructed for each photometric band. In the I-band, a 



M. Reyniers et al.: The rotation and coma profiles of comet C/2004 Q2 (Machholz) 



clear m--\ for p<10 4 km and m = -1 .5 for p>10 4 km was 
found, in agreement with a steady state coma without creation 
or destruction of particles. From this regime change, a termi- 
nal ejection velocity of the grains v gr = 1609 m s was derived. 
In the other filters, smaller m-values were found, due to the con- 
tamination of strong gas features in those bands, while the I filter 
is relatively free from those emissions. 
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